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Abstract 

If an ultrarelativistic charged particle channels inside a single crystal with periodically bent crystallographic planes, it emits hard 
electromagnetic radiation of the undulator type. Due to similarity of its physical principles to the ordinary (magnetic) undulator, 
such a device is termed as the crystalline undulator. Recent development of a new Monte Carlo code ChaS made possible a detailed 
simulation of particle channeling and radiation emission in periodically bent crystals. According to recent findings, energy of the 
electron beam below 1 GeV is sufficient to observe the undulator effect in a periodically bent crystal. Even more exciting results 
were obtained for a crystalline undulator whose bending period is shorter than the period of the channeling oscillations and the 
bending amplitude is smaller than the width of the planar channel. Such a crystalline undulator is far superior to what was proposed 
previously. It allows for a large effective number of undulator periods. Therefore, it is predicted to emit intense undulator radiation 
in the forward direction. A narrow undulator peak is seen for both positron and electron beams. Using positrons is, however, more 
desirable because in this case the intensity of the undulator radiation is higher while the background is lower. 

Keywords: channeling, radiation by moving charges, crystalline undulator, channeling radiation, electron and positron beams, 
synchrotron radiation sources, Monte Carlo method, strained layer superlattice 
PACS: 61.85.+p,41.60.-m,41.75.Fr,02.70.Uu 
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The idea of the crystalline undulator (CU) was proposed 
more than three decades ago It was suggested 

to use a periodically bent crystal to generate hard electromag- 
netic radiation of undulator type. Ultrarelativistic electrons or 
positrons were supposed to channel through such a crystal fol- 
lowing the sinusoidal shape of the bent crystallographic planes. 
Due to nearly harmonic transverse oscillations of the parti- 
cles, the electromagnetic waves emitted in the forward direction 
were expected to have a narrow spectral distribution, similarly 
to an ordinary (magnetic) undulator (5J y, [J . 

The extremely strong electromagnetic fields inside the crys- 
tal can steer the beam particles much more effectively than it 
would be possible even with the best superconductive magnets. 
Therefore, the period of the crystalline undulator can be made 
orders of magnitude smaller than that of the magnetic undula- 
tor. Hence, the crystalline undulator is potentially capable of 
generating hard X rays and gamma rays. 

There is, however, a price to pay for taking advantage of the 
crystalline field. In contrast to the ordinary undulator, the parti- 
cles move in a dense medium instead of vacuum. They experi- 
ence random collisions with crystal constituents and, therefore, 
they emit bremsstrahlung which may constitute a substantial 
background. Coherent effects contaminate the spectrum even 
stronger. In addition to undulator oscillations, the particle has 
to perform channeling oscillations around the minimum of the 
planar potential (similarly to its motion in a straight crystal). 
As a result, channeling radiation is present in the spectrum of 



the crystalline undulator in addition to undulator radiation and 
bremsstrahlung. 

Channeling radiation has a lot in common with undulator ra- 
diation Jit]. Therefore, even a straight crystal can be used to 
build a source of hard photons. However, this approach has a 
disadvantage. Because the shape of the transverse potential is 
not parabolic, the transverse motion of the particles is not har- 
monic. As a result, the spectrum of the channeling radiation is 
broader than that of the undulator radiation, especially in the 
case of negatively charged projectiles^ 

An interesting phenomenon of narrowing the planar channel- 
ing radiation peak has been observed for positrons H 1 3l 1 1 411 . The 
reason is a partial compensation of the potential inharmonic- 
ity by deviations from the dipole approximation in the radia- 
tion emission 1 15]. This takes, however, place only at a certain 
'magic' beam energy. This means that the position of the nar- 
row peak is fixed for every crystal channel, i.e. the radiation 
frequency can be varied only within a very limited interval. 
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1 Another difference between channeling radiation and undulator radiation 
is that the latter one can be coherent if the particle beam is modulated (bunched) 
in the longitudinal direction with the period equal to the wavelength of the un- 
dulator peak. In this case the different particles radiate electromagnetic waves 
with nearly the same phase. Therefore, the intensity of the radiation becomes 
proportional to the number of particles squared (in contrast to a linear propor- 
tionality for an unmodulated beam). This way of producing coherent radiation 
is utilized in free electron lasers [9] (see e.g. [ 10] for a modern review). Simi- 
lar effect should be observed in a crystalline undulator if it is fed by a bunched 
particle beam flUllal . In contrast, channeling radiation will not become coher- 
ent even in the case of a modulated beam because the phases of the channeling 
oscillations are arbitrary and random. 
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The CU approach does not have such restrictions. The un- 
dulator radiation peak can be shifted within a wide range of 
values by varying the beam energy and the bending period. But 
one has to put up with the contamination of the spectrum by 
channeling radiation unless one finds a way to get rid of it. 

The radiation background is not the only challenge faced by 
the CU. It was realized at the very beginning [4] that the ef- 
fective length of the crystalline undulator is limited by the at- 
tenuation of photons in the crystal medium. This limitation is, 
however, essential if rather soft photons, ha> < 100 keV are to 
be produced. 

A more severe restriction on the effective length of the bent 
crystal is imposed by the dechanneling phenomenon [16]. The 
incoherent bremsstrahlung is not the only undesirable effect of 
random scattering of the projectiles by the crystal constituents. 
Due to the random collisions, the transverse energy of the chan- 
neling particle fluctuates. Positive fluctuations are, however, 
more likely than negative ones. Therefore, the projectile gains 
on average the transverse energy. If the latter exceeds the 
height of the interchannel potential barrier, the particle leaves 
the channel 11711 . Starting from this point, it does not follow the 
shape of the channel and, consequently, it does not emit undu- 
lator radiation. For this reason, the effective number of undula- 
tor periods N u is limited by the average length Ld at which the 
dechanneling takes place: 



N ~^ 

u V 



(i) 



where A u is the bending period of the crystalline undulator. The 
number of undulator periods has to be large, N u » 1, to ensure 
a narrow spectral distribution of the undulator radiation. There- 
fore, the following condition has to be satisfied by the undulator 
period 

« U (2) 

In its initial form, the idea of the CU was based on the as- 
sumption that the projectile should follow the sinusoidal shape 
of the bent crystallographic planes and perform, at the same 
time, channeling oscillations around the central plane of the 
channel. This implied that the bending period of the undulator 
had to be much larger than the period of channeling oscillations 



A u » Ac 



(3) 



In addition, the undulator bending amplitude a u has to be 
much larger than the typical amplitude of channeling oscilla- 
tions a c . This condition can to be rewritten in the form 1118 1 



a a » d 



(4) 



due to the fact that a c < d/2, where d is the channel width 
(the distance between the bent crystallographic planes that form 
the channel). The strong inequality © has to be satisfied to 
ensure a higher intensity of the undulator radiation relative to 
the channeling radiation. 

In the following, the crystalline undulator satisfying condi- 
tions (0 and (H]i will referred to as LALP CU (large amplitude 
and long period crystalline undulator). The complete list of 



conditions that have to be satisfied by the parameters of LALP 
CU can be found in |[l9l]0 

One more condition is relevant to the present discussion. It 
ensures a stable channeling of the projectile in the periodically 
bent crystal of LALP CU ||2>l3l- Initially, similar condition was 
obtained for channels with constant curvature ll20ll . It is conve- 
nient to write it down in the formal 

F c f . t a„E 



\>C = 



UL 



= 4n 2 



(5) 



Here C is the centrifugal parameter IU9II . F c f is the maximal 
centrifugal force acting on the projectile in a the periodically 
bent channel, U' milx is the maximal force that keeps the particle 
in the channel and E is the energy of the projectile. 

Conditions (0, and (0 are difficult to satisfy simulta- 
neously. In particular, they cannot be satisfied in the case of 
electron beam of moderate energy, E < 1 GeV. For instance, 
if one sticks with (0, the most favorable conditions for observ- 
ing undulator radiation from E = 855 MeV electrons in the 
Si(110) channel are A u » 4 fim and C ss 0.3. This corresponds 
to a a ~ 0.84 A (a u /d ~ 0.44), i.e. condition © is broken. Tak- 
ing into account that the dechanneling length L d a 8.3 fim even 
for the straight Si(110) channel [21]Q one sees that the strong 
inequality (0 is not satisfied either. As a result, the undulator 
peak is rather small and not very sharp while the spectrum is 
dominated by the channeling radiation (see figure0. 

Still, it is remarkable that the undulator effect is predicted 
to be detectable even using electron beams of moderate ener- 
gies. This awakes expectation that a successful proof of princi- 
ple experiment can be done in the nearest future, for example, 
at Mainz Micro trone (MAMI). 

The LALP CU conditions can be fulfilled for positron beams 
(see 11911 and references therein) and for high energy, E > 10 
GeV, electrons [23]. 

An example of a spectrum of a positron-based LALP CU is 
shown in figure [2] The undulator peak is narrow and is sub- 
stantially higher than the channeling radiation maximum. Nev- 
ertheless, the total energy of channeling radiation (integrated 
over the photon energy interval 0.8 MeV < fuo < 1.2 MeV) 
exceeds that of undulator radiation. This is even more true for 
high energy electrons (see e.g. figure 4 of 123J). Moreover, 
the channeling photons are harder than the undulator ones and, 
therefore, they cannot be easily screened out. This may cause 
serious problems for many potential applications. 



2 It was argued in 1T9I1 that condition {4} is necessary to ensure separation of 
the undulator radiation and the channeling radiation in the spectrum. In fact, it 
is not so. The separation is ensured by y}. The latter inequality does follow 
from l|4j when it is supplemented by {5], but (3) can be fulfilled even if a u < d. 
However, the intensity of the undulator radiation will be small in this case. This 
is the reason why l|4j is necessary in the case of LALP CU. 

3 The centrifugal parameter can be also expressed in the form C = R c /Rmm, 
where R c = £/J/^ ax in the critical radius of the channel 1 20] (also known as 
Tsyganov's radius) and R m j n = A^I(An 2 a a ) is the minimal curvature radius of a 
sinusoid with the period A a and the amplitude a u . 

4 Dechanneling is defined in 12111 as crossing the channel boundary. Other 
authors (see e.g. 12211 ) define dechanneling as rising the energy of the channel- 
ing oscillations above the interchannel barrier. The value of the dechanneling 
length depends only slightly on the definition that is used for its computation. 
The difference does not exceed a few hundreds of nanometers. 



2 




Figure 1 : Simulated spectrum of the radiation emitted by E = 855 MeV elec- 
trons channeling in a 24 /im long LALP CU with the bending period A u = 4 fjm 
and the centrifugal parameter C = 0.3. The undulator radiation peak is seen at 
hw s* 1.6 MeV. The spectrum is dominated by channeling radiation (the broad 
maximum at 5 MeV < Imj < 9 MeV ). 



The results of Monte Carlo simulations shown in figures 
Q] and |2 were obtained with a new computer code ChaS 
(Channeling Simulator). The code performs a 3D simulation 
of particle motion in the crystal and calculates the spectral and 
angular distribution of the emitted radiation. 

The previously existing channeling codes can be divided into 
two main categories. There are codes based on the continuous 
potential approximation J3 |H S |27l H H H and those 
considering binar y co llisions of the projectile with the crystal 
atoms i3Ul32ll33ll34il35ll . In the latter case, the atom is usually 
taken as as a whole, i.e. incoherent collisions with atomic elec- 
trons are ignored. In contrast, the algorithm of ChaS is based 
on the binary collisions of the projectile with target electrons 
and target nuclei. This novel approach is especially beneficial 
in the case of negatively charged projectiles, that have to cross 
the crystal plane during the channeling process. Indeed, the 
continuous potential approximation becomes inaccurate in the 
vicinity of the atomic nuclei while the electron density near the 
crystal plane is much higher than the average one. 

The previous version of the code used the electron distribu- 
tion calculated from Moliere's parametrization of the atomic 
potential. The obtained results were published in 112 ll 13611 . 
They demonstrated reasonable agreement with experimental 
data 112 ill . The present version of ChaS employs a more ef- 
ficient and robust algorithm for the calculation of the emitted 
radiation. In addition, it has an option of using the first princi- 
ple distribution of electrons in the crystal instead of Moliere's 
parametrization. The first principle distribution is calculated 
within the density functional theory using the computer code 
ABINIT f37tl. 



Figure 2: Simulated spectrum of the radiation emitted by E = 500 MeV 
positrons channeling in a 350 ^m long LALP CU with the bending period 
A u = 23.3 fim and the centrifugal parameter C ra 0.12 (^, = 19.2 A, a u /d as 10). 
The undulator radiation peak is seen at hoi a 88 keV. The parameters of the 
crystalline undulator are the same as in right bottom panel of figure 8 of (Hi. 
In contrast to 1131 . the present simulations take into account the incoherent 
background. 



The code ChaS in its present version ist most suitable for 
the analysis of the channeling of electrons and positrons with 
energy E in the range from a few hundreds of MeV to a few 
GeV with the emition of not very soft: ha> > 0.5 MeV and at 
the same time not too hard hat <k E photons. The algorithm 
takes into account all the physics that is relevant to this domain. 
A number phenomena are neglected in the model (e.g. influence 
of the crystal medium on the emission and propagation of the 
radiation, quantum effects in the motion of the projectile, losses 
of the projectile energy due to emission of photons, a shift of 
the photon energy due to recoil etc.) They are expected to be 
small and do not influence the results substantially J21 , 381. 

For simplicity, the emittance of the particle beam was ne- 
glected in all simulations presented in this article. The particles 
were assumed to enter the crystal at zero angle to the crystallo- 
graphic planes. This is a reasonable approximation in the case 
of channeling experiments with high quality electron beams 
1 39]. It may not be the case for position beams, but the beam 
divergence depends on details of the experimental conditions 
whose analysis is out of the scope of the present contribution. 



It was suggested recently [40] that conditions and (0 are 
not necessary. In fact, an intense source of hard photons with 
a narrow spectral distribution can can be created if both condi- 
tions are violated. 

First, let us reanalyze the reasons behind condition (|4). It is 
needed to make sure that the spectrum is dominated by the un- 
dulator radiation rather than by the channeling one. However, 
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the amplitude of undulator oscillations has to be much larger 
than that of the channeling oscillations only if the frequency 
of the undulator radiation w u is smaller than the frequency of 
channeling oscillations <y c . Indeed, the energy radiated in a cer- 
tain direction (the forward direction in the present case) by a 
moving particle in the dipole approximation has the following 
dependence on the transverse oscillation amplitude a and the 
radiation frequency a>: 



dS 



da>dQ. 



a 2 a> 4 . 



(6) 



Here dQ. is the differential of the solid angle and 9 is the angle 
between the direction of the radiation emission and the average 
direction of the particle motionjf] Therefore, condition © is 
not necessary, i.e. the amplitude of the undulator bending can 
be smaller than the channel width, 



o u < d, 



(7) 



if the frequency of the undulator radiation is considerably larger 
than that of the channeling radiation 



(8) 



To fulfill this condition, the period of the crystal bending A u 
has to be much smaller than the smallest period of channeling 
oscillations Ac: 

A a « Ac. (9) 

The last inequality violates condition (0. This can be seen from 
the following consideration. The period of the channeling os- 
cillations can be estimated by 



2n 



u"(oy 



(10) 



where t/"(0) is the second derivative of the transverse potential 
energy with respect to y in the point of its minimum y — (the 
axis y is perpendicular to the channel boundaries). Taking into 
account thao 

U' max <U"(0)d (11) 
in combination with ( ITOb and (O one obtains from (|5j 



1 > C » 



(12) 



The bending amplitude of the crystalline undulator a u cannot 
be much smaller than the channel width d otherwise it becomes 
comparable to (or even smaller than) the amplitude of thermal 
vibrations of the atoms in the crystal. Clearly, the undulator 



5 The dependence {6) becomes obvious from the proportionality of the ra- 
diated energy to the particle acceleration squared. For transverse harmonic 
oscillations, the transverse acceleration is proportional to the oscillation fre- 
quency squared and linearly proportional to the oscillation amplitude. Hence 
the dependence {6) is obtained. 

Expressions j 10) and lilt would be exact equalities in the case of parabolic 
potential. For a real potential, the second derivative is not constant and the 
maximum value of the force is reached at \y\ < d, hence )10t is an approximate 
equality and )1U is an inequality. 
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Figure 3: Simulated trajectories of a positron (upper panel) and an electron 
(lower panel) with energy E = 855 MeV channeling in a crystalline undulator 
with a small amplitude, a u = 0.4 A, and a short period, A u = 400 nm. The 
projectile does not follow the shape of the bent crystallographic planes (the 
thick wavy lines). It performs channeling oscillations with roughly the same 
period as in a straight crystal. The effect of crystal bending on the shape of 
tr aject ories is barely seen. The figure is a modified version of Fig. 1 published 



effect will not be present in this case. If a u is comparable to d, 
the two inequalities of (fT2l become incompatible. One might 
expect that it destroyes the undulator effect, but, fortunately, it 
does not. Condition (0 is, in fact, not applicable in the case of 
SASP CU (small amplitude © and short period (0 crystalline 
undulator). 

In figure [3] simulated trajectories of a positron and an elec- 
tron channeling in a SASP CU are plotted. One sees from the 
figure that the particles do not follow the shape of the bent crys- 
tallographic planes. Therefore, formula (|5} is irrelevant. The 
channeling process is still present. The particle motion can 
be roughly considered as if it were governed by a continuous 
potential averaged over the oscillations of the plane. In other 
words, it is similar to the channeling in a straight crystal but the 
average transverse potential is somewhat different. 

Still, a more thorough consideration reveals that the particle 
performs also transverse oscillations with the undulator period 
A u . The amplitude of these oscillations a is much smaller than 
the bending amplitude a u . Therefore, it is difficult to see the 
modification of the trajectories induced by the crystal bending 
in figure [3] Nonetheless, the corresponding Fourier harmon- 
ics are present and they are seen in the spectra of the emitted 
radiation shown in figures |4] and [5] 

The undulator radiation peak is higher than the channeling 
radiation maximum despite of the amplitude a of the undulator 
oscillations of the projectile being much smaller than the am- 
plitude of its channeling oscillations a c . This seemingly para- 
doxical fact can be easily understood if relations © and © are 
taken into consideration. 

The undulator peaks are narrow and well separated from the 
channeling radiation. The undulator and channeling radiation 
peaks have approximately the same absolute width. Nonethe- 
less, the relative width of undulator peaks is much smaller due 
to®. 

Due to small A u , condition ([TJ can be satisfied for SASP CU 
even if it is used with a moderate energy electron beam. This 
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Figure 4: Simulated spectra of the radiation emitted by E = 855 MeV positrons 
channeling in a 12 j/m long crystalline undulator with a small bending ampli- 
tude for four short periods. The spectra are obtained by averaging over about 
1000 simulated trajectories. The undulator radiation peaks are higher and are 
centered at a much larger photon energy than the corresponding channeling 
radiation (ChR) peaks. 



is a big practical advantage since electron beams are usually of 
higher quality and are less expensive than positron ones. Still, 
positrons are more preferable. As one sees from figures 0and |5] 
they allow for a higher intensity of the undulator peak accom- 
panied by a lower background. 

The undulator radiation of S ASP CU is much harder than the 
channeling radiation (8j. It is an important advantage of SASP 
CU over LALP CU. It is usually much easier to get rid of a 
soft photon background and preserve the hard part of the spec- 
trum than to do the opposite. For example, a detector can be 
made sensitive to hard photons, but be screened from (or made 
insensitive to) soft photons. It is difficult and sometimes even 
impossible to do vice versa. For this reason, SASP CU is ex- 
pected to be much more suitable for many potential applications 
than LALP CU. 

Due to its much smaller bending period, SASP CU can pro- 
duce by about two orders of magnitude harder photons when 
used with a beam of the same energy as LALP CU. Or, in other 
words, SASP CU will require a much smaller and, therefore, a 
much less expensive accelerator than the one which would be 
needed for the production of radiation of the same frequency 
with LALP CU. 

Hence, the crystalline undulator that violates conditions © 
and (0 has a number of advantages with respect to LALP CU. 

From the technological point of view, SASP CU is more chal- 
lenging than LALP CU. There exist at least four technologies 
suitable for the fabrication of LALP CU. The oldest idea of us- 
ing ultrasonic waves |2j, |4J] is, unfortunately, still waiting for 
its experimental implementation. Two other technologies uti- 
lize the idea of imposing periodic stresses on the surface of the 



Figure 5: The same as in figure |4]but for electrons. The undulator radiation 
peaks are lower than in figure|3]by a factor of ~ 3 and the background is higher. 



crystal sample. These are making regularly spaced groo ves on 
the crystal surface either by a diamond blade ll4ll l42ll or by 

means of laser-ablation [43] and deposition of periodic S13N4 

i — I 

layers onto the surface of a Si crystal [42]. Finally, there is a 
way to create periodically varying stresses in the crystal vol- 
ume by growing a crystal with periodically varying chemical 
composition. The most mature technology is the creation of 
strained layer [44] superlattices by periodically varying germa- 
nium concentration x in a Sii-^Ge* crystal with a periodically 
varying Ge content x l45l l46ll . 

The potential applicability of ultrasound in the case of SASP 
CU requires further investigations. The methods based on sur- 
face stresses cannot be applied because they require the trans- 
verse dimension of the crystal to be of the order of the bend- 
ing period B47I1 . The latter is smaller than 1 fim in the case 
of SASP CU. Only the last approach, the growing of Sii_ A Ge A 
strained layer superlattices, is suitable for the fabrication of 
SASP CU. This technology has been already used for manufac- 
turing the LALP CU that is being used in ongoing experiments 
at Mainz Micro tron ll48ll . It was demonstrated recently [49] 
that a Sii_ A Ge A crystalline heterostructure can be grown with 
the parameters that have been used in the simulations presented 
in figures[4]and[5] Moreover, the strained layer crystal with pa- 
rameters of SASP CU was predicted to be stable agains misfit 
dislocations, in contrast to LALP CU which is only metastable. 

In conclusion, the crystalline undulator with a small ampli- 
tude and a short period can be created and it is predicted to be 
far superior with respect to LALP CU. 

In the present contribution, production of undulator radia- 
tion with photon energy in the range of tens of megaelectron- 
volts is considered. These results are important because the 
hard photon range is unattainable for the present state-of-the-art 
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synchrotron radiation sources. It would be, however, interest- 
ing to consider production of softer photons, in the range of a 
few hundreds or even tens of kiloelectronvolts. This domain is 
on the edge of the capability of the presently existing and con- 
structed facilities ifToL [50h - This facilities are, however, unique 
and very expensive. Due to the fact that SASP CU requires 
a much smaller accelerator than a conventional undulator for 
production of photons of the same energy, it has a potential to 
offer a much less expensive solution. If it were possible to use 
SASP CU for production of soft X rays, such devices could be 
made affordable even to medium size university labs or hospi- 
tals. A theoretical investigation of the low energy SASP CU is, 
however, more challenging, because quantum effects definitely 
cannot be neglected in this case. 

Even more exciting future task of this field of science is ex- 
ploring the possibility to produce coherent radiation with SASP 

cuGj. 
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